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Full QCD with dynamical Wilson fermions on a 24 3 x 40-lattice - 
a feasibility study* 
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The investigation of light sea-quark effects in lattice QCD with dynamical Wilson fermions requires both larger 
physical volumes and finer lattice resolutions than achieved previously. As high-end supercomputers like the 
512-node APE Tower provide the compute power to perform a major step towards the chiral limit (T^L), we have 
launched a feasibility study on a 24 3 x 40 lattice. We approach the chiral limit-while refining the resolution-, 
using the standard Wilson fermion action. Following previous work, our Hybrid Monte Carlo simulation runs at 
p — 5.6 and two K-values, 0.1575 and 0.158. From our study, we are confident that, for the APE Tower, a realistic 
working point has been found corresponding to a volume of 2 fm 3 , with chirality characterized by ~ 5.6. 



1. INTRODUCTION 

Simulations of full QCD with Wilson fermions 
at zero temperature so far have been carried out 
on lattices of size < 16 3 x 32, physical volumes < 
1.5 (fm) 3 and ratios of ^ > 0.6 §. The latter 
quantity is a monitor for the closeness to the chi- 
ral point. It has been demonstrated that a sta- 
tistically significant full QCD (reference) sample 
can be generated in one year's runtime on a 256- 
node APE computer, at ^ = 0.71. This however 
corresponds still to a rather heavy quark mass 
m q : by use of chiral perturbation theory we can 
mock a fictitious "physical" pseudoscalar meson 
containing two strange quarks, with mass ratio 

of the size quoted, ^ m V?™5l = o.69. Thus, 
in order to quantify light sea quark effects in full 
QCD, one would rather prefer to work on larger 
volumes that accommodate a large 7r-corrclation 
length both in physical and lattice units. This 
clearly asks for simulations on lattices > 16 3 and 
f3 > 5.6. 

In this note, we describe the T%L project, 
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which is geared to push QCD simulations with 
standard Wilson fermions further towards the chi- 
ral limit, i.e. beyond < 0.6 and at appropriate 
volumes. 

The 512-node APE Tower offers sufficient 
memory to handle a 24 3 x 40 lattice. With its 
CPU-power it can drive an optimized HMC at 
sufficient speed (?) to increase the lattice size by 
more than a factor of 4 compared to the previ- 
ous standards described above, (ii) to go more 
chiral, i.e., cope with worse conditioned fermion 
matrices. 

In our exploratory study, we are guided by the 
experiences described in taking advantage 

of algorithmic achievements such as improved in- 
verters J4| and new parallel preconditioning tech- 
niques g. We shall report on a Hybrid Monte 
Carlo simulation on a 24 3 x 40 lattice at (3 = 5.6 
and two K-values, 0.1575 and 0.158. 

2. HOW CHIRAL? 

The 24 3 x 40 lattice allows to increase ^ by a 
factor of 1.5 compared to Ref. J(|, which should 
suffice to target for ^ in the range of .5 - .6. 
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Concerning physical volumes and scales we ben- 
efit from the increasing A/3-shift H as we go to 
smaller bare quark masses and choose (3 = 5.6. 
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Figure 1. Fixing Ksea and estimate for k c 
compared to the results from Ref H. 



For the determination of n seai we extrapolated 
the relation m q a = | (j^ — on the data set of 

Ref. § to m q a = 0.023, cf. Fig. ga, where is 
estimated from the mass trajectory to be about 
.23 of the spatial extension. This value is small 
enough to protect us from finite size effects. We 
thus will work at k = 0.1580. 

To put this parameter choice into perspective, 
we sketched in Fig. |l|b the approach of the critical 
« at fixed K SO a, tc,v; Hj towards k c , the locus of 
which is ftgea = n c (the diagonal line). The cross 
marks the current estimate for our working point. 
Notice that our parameter choice, « S ca = 0.158, 
appears to be reasonably positioned within the 
'chirality gap'. 

In the lattice discretization of the fermionic ac- 
tion, we switched from the usual o/e representa- 
tion to the full representation of M = 1 — kD, 
employing a new SSOR preconditioning scheme 
HgD- In particular on the APE machine, this 
method offers an overall gain of 100 % in execu- 
tion time, as seen from Tab. g. 



Algo 




K = 0.1575 re = 0.1580 


o/e 


t/s 


8200 


SSOR 


t/s 


3800 9100 



Table 1 

Average time to generate 1 trajectory on the 



APE100 Tower. 



3. SIMULATION 



We have tuned the HMC timestep to achieve 
acceptance rates larger than 60 %. For the SSOR 
scheme with twice as many degrees of freedom 
as in the o/e case, we chose T = 0.5. For this 
trajectory length, in the production runs, the 32 
bit machine precision induces a reversibility error 
6(AS) in the range of just 2 % of the average 
AS for an inversion residue of r = 10~ 8 . This 
is due to local computations, global summations 
being carried out in emulated double precision 
arithmetic. It should be said that the impact of 
this error of AS onto the canonical distribution 
deserves further attention. The chosen HMC run 
parameters are given in Tab. ^. 



Algo 


K 


T 


dt 


acc/% 


r 


o/e 


0.1575 


1 


0.008 


70 


10 _s 


SSOR 


0.1575 


0.5 


0.004 


72 


lO" 8 


SSOR 


0.1580 


0.5 


0.004 


66 


io- s 



Table 2 

Parameters of the HMC simulation. 



During the thermalization phase we carefully 
approached the lowest quark mass in a near adia- 
batic fashion, to protect the system from oscillat- 
ing through the shielding transition. We forked 
the run into two K-branches after an initial ther- 
malization of 480 trajectories. The production 
status reported is given by 830 trajectories (out 
of which 350 are thermalized) for k = 0.1575 and 
1500 trajectories (out of which 750 are thermal- 
ized) for k = 0.1580. 
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4. PRELIMINARY RESULTS 

We have performed first tentative measure- 
ments of the autocorrelation C(t) of the plaque- 
tte. In Fig. pi we plot the autocorrelation func- 
tion for the plaquette at re = 0.158. On the large 
lattice volume we can profit from a substantial 
self averaging effect suppressing the fluctuations 
of the plaquette as well as of other intensive quan- 
tities. The plot indicates that the autocorrela- 
tion times come out surprisingly small and might 
settle well below r m + = 50 for the plaquette. 
We mention that the autocorrelation function for 
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Figure 2. Autocorrelation function of plaquette. 



light meson masses looks similar. 

Using 22 configurations drawn from a sample of 
600 trajectories we have computed the potential 
following Ref. S. We performed a 30 step APE 
smearing and evaluated the potential for a time 
extension of 5 where a plateau in the local mass 
is emerging. We have fitted for a " string tension" 
in the range up to 1 fm. In Tab. |^ we quote a 
preliminary estimate for the ensueing scale and 
physical lattice volume V s . 

We monitored local meson masses to position 
the run with respect to chirality: on a sample of 
19 configurations we retrieve a rough first guessti- 
mate of m^a and m p a, see Tab. |3J. Our findings 
suggest ^ to be 0.56(4). 



or 1 V s m-x SL ra p SL ^ 

2.37(1) GcV 2 fm 0.178(5) 0.32(2) 0.56(4) 



Table 3 

Results for the lattice scale from potential, and 
masses in lattice units. 



5. CONCLUSIONS AND OUTLOOK 

In the T^L-feasibility study, we find that ^f— 

ratio appears to reach the target region indeed, 

f i 
where — = 5.6 < 0.25 x V s 3 . The lattice rcsolu- 

tion is increased (with respect to the small lattice 

results) to a -1 = 2.37 GeV. We are encouraged 

by the observed autocorelation times and expect 

> 50 independent configurations from 8 months 

future runtime on APE Towers. 
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